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; Abstract 

Q ■ Using the bond-propagation algorithm, we study the finite-size behavior of the critical two- 

^ ■ dimensional Ising model on a finite triangular lattice with free boundaries in five shapes: triangle, 

i rhombus, trapezoid, hexagon and rectangle. The critical free energy, internal energy and specific 

■ heat are calculated. The accuracy of the free energy reaches 10~^^. Based on accurate data on 
> I 

g ■ several finite systems with linear size up to = 2000, we extract the bulk, surface and corner 

. parts of the free energy, internal energy and specific heat accurately. We confirm the conformal 

Q \ field theory prediction of the corner free energy to be universal and find logarithmic corrections 
in higher order terms in the critical free energy for the rhombus, trapezoid, and hexagon shaped 

■ systems, which are absent for the triangle and rectangle shaped systems. The logarithmic edge 
, corrections due to edges parallel or perpendicular to the bond directions in the internal energy are 

o- 

PsJ ■ found to be identical, while the logarithmic edge corrections due to corresponding edges in the free 

(N. 



X: 



energy and the specific heat are different. The corner internal energy and corner specific heat for 
angles vr/S, tt/2 and 27r/3 are obtained, as well as higher order corrections. Comparing with the 
corner internal energy and corner specific heat previously found on a rectangle of the square lattice 
(Phys. Rev. E. 86 041149 (2012)), we conclude that the corner internal energy and corner specific 
heat for the rectangle shape are not universal. 
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I. INTRODUCTION 



The finite-size scaling theory, introduced by Fisher, finds extensive apphcations in the 
analysis of experimental, Monte Carlo, and transfer-matrix data, as well as in recent theoret- 
ical developments related to conformal invariance [l-4|. It becomes of practical interest due 
to the recent progresses in fine processing technologies, which has enabled the fabrication of 



nanoscale materials with novel shapes 



-|7|]. Exact solutions have been 



playing a key role 



in determining the form of finite-size scaling. Ferdinand and Fisher js] pioneered on the 
two-dimensional (2D) Ising model on a finite size square lattice, which extended Onsager's 
exact solution jsj and stimulated the ideas of finite-size scaling. Since then, exact results of 
the model on finite size lattices with various boundaries have been studied intensively 



Q. Detai 
condition 



ed knowledge has been obtained for the torus case [l2|, for the helical boundary 



13| . for the Brascamp-Kunz boundary condition |l^.]l5[| and for an infinitely long 



cylinder |l6| . The exact solution on the triangular lattice has also been studied intensively 

HQ. 

However for the 2D Ising model the exact solution with free boundaries, i.e. free edges and 
sharp corners, is still missed. As we know, the Bethe ansatz is a powerful technique to solve 
the 2D Ising model, but this technique would not allow to place the system on a rectangle 
with free boundary conditions on top and bottom, since the expression of the boundary 
state in terms of the Bethe eigenvectors is unknown, which is a famous unsolved problem 
in statistical mechanics. Although there are Monte Carlo and transfer matrix studies on 
this problem 19|, |20|, the accuracy or the system sizes of the results are not enough for 
extracting the finite size corrections. Meanwhile, for 2D critical systems, a huge amount of 
knowledge has been obtained by the application of the powerful techniques of integrability 



and conformal field theory (CFT) 



21 



-|23|. Cardy and Peschel predicted that the next 



subdominant contribution to the free energy on a square comes from the corners 



3i , which is 



universal, and related to the central charge c in the continuum limit. Kleban and Vassileva 



2l| extended the result to a rectangle. These results are consistent with the conjectured 



exact analytic formula for the Ising model on a square lattice 24l |. 

Several years ago an efficient bond propagation (BP) algorithm was developed for com- 
putin g th e partition function of the Ising model with free edges and corners in two dimen- 



sions 



25 



26| . Making use of this algorithm, we recently determined numerically the exact 



partition function of the Ising model on the square lattice with rectangle shape and free 
boundaries {2^. We not only confirmed the CFT predictions, but also found logarithmic 
corrections due to corners in the internal energy and specific heat. 

In present paper we apply the BP algorithm to study the Ising model on finite triangular 
lattices with free boundaries in five different shapes, focusing on how the shape affects the 
finite-size scaling. The five shapes are triangle, rhombus, trapezoid, hexagon and rectangle, 
as shown in Fig. [1] Based on accurate data on a sequence of finite systems with linear size up 
to = 2000, we extract the bulk, surface and corner parts of the free energy, internal energy 
and specific heat accurately. We verify the conformal field theory prediction of the corner 
free energy and find logarithmic corrections in higher order terms in the critical free energy 
for the rhombus, trapezoid, and hexagon shaped systems, which are absent for the triangle 
and rectangle shaped systems. The logarithmic edge corrections due to edges parallel or 
perpendicular to the bond direction in the internal energy are found to be identical, while 
the logarithmic edge corrections due to the corresponding edges in the specific heat are 
different. The corner internal energy and corner specific heat for angles 7r/3, 7r/2 and 27r/3 
are obtained, as well as higher order corrections. Comparing with the previous found corner 
internal energy and corner specific heat on a rectangle of the square lattice {2^ , we conclude 
that the corner internal energy and corner specific heat for the rectangle shape are not 
universal. 

Our paper is organized as follows: In section [Tll we briefly describe the BP algorithm 
used here. We present in sec. [Tlll our main results and discussion. We conclude in sec. HVl 

II. METHOD 

The partition function of the Ising model on the 2D triangular lattice is 



where the nearest neighbor couplings are dimensionless and /3 is the inverse temperature. 
This partition function for finite triangular lattice with five different shapes and open bound- 



0.274653072167 ■■■ . Figure [T] shows the five shapes: triangle, rhombus, trapezoid, hexagon 
and rectangle. The linear size of a finite lattice is defined as the length of edges in the tri- 
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aries is calculated with the BP algorithm 25|] at the exact critical point Pc 



I ln(3) 
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FIG. 1. (a) The triangle-shaped triangular lattice with = 4. (b) The rhombus-shaped lattice 
with = 4. (c) The trapezoid-shaped lattice with = 3. (d) The hexagon-shaped lattice with 
= 3. (e) The rectangle-shaped lattice with A^ = 5. Three bond directions and the perpendicular 
direction are shown (see text). 



angle, rhombus and hexagon cases, of which the length of edges are equal. For the trapezoid 
shape, N is the shortest edge, of which the lengths of the three short edges are required to 
be equal. For the rectangle, is defined as the length of the bottom edge, and the number 
of layers is required to be A^, which means the actual geometrical vertical length is N\/3/2. 
According to the finite-size scaling the system size should be the actual geometrical 
length, not the number of layers. Therefore the aspect ratio of the rectangle, we consider 
here, is \/3/2 rather than 1. 



The B 



algorithm for the Ising model on the triangular lattice has been described in 



detail in 26|]. For the triangle-shaped and rhombus- shaped lattice, this algorithm can be 



applied directly. However, for the trapezoid-shaped, hexagon-shaped and rectangle-shaped 
lattice, the inverse of the BP series reduction should be introduced, which corresponds to 
generating a new spin between two spins such that 



J120-1 (72 



,5F+Jao{<Ji+a2) 



as illustrated in Fig. |2]^a). It is convenient to use variables j = e 



,Jl2 



(2) 

e--^i2 and 5f = e^^ . 



Then we get the solution 6f = 712/2 and j = \ ^ — \ — 1. 

In each step of the BP algorithm the transformation is exact. The numerical accuracy is 
only limited by the machine's precision, which is the round-off error 10~^^ in the quadruple 




FIG. 2. (a) Inverse of the BP series reduction. See the formulation for this transformation in text. 
(b)-(g) are schematic of the algorithm for the trapezoid-shaped lattice. From (b) to (c), the inverse 
of BP series reduction is applied to the diagonal bond ( thick line) at the bottom of the trapezoid. 
From (c) to (d), this operation is applied to another diagonal bond (thick line). From (d) to (e) 
bond propagation operations are applied. From (e) to (f), the inverse of BP series reduction is 
applied again, (f) to (g) Those diagonal bonds are eliminated by the usual BP procedure. 

precision. The BP algorithm needs about A^^ steps to calculate the free energy of a. N x N 
lattice (much faster than other numerical method). Therefore the total error is approxi- 
mately N^/"^ X 10^'^^. This estimation has been verified in the following way: We compared 
the results obtained using double precision, in which there are 16 effective decimal digits, 
and those using quadruple precision. Because the latter results are much more accurate than 
the formal, we can estimate the error in double precision results by taking the quadruple 
results as the exact results. We thus found that the error is about A^3/2 X 10-16. In our 
calculation, the largest size reached is = 2000, the round-off error is less than 10^^^. 

The free energy density, internal energy per spin and specific heat density are calculated 
at the critical point according to 



InZ _dl _ 2^ 



(3) 



respectively, where S is the number of spins on the lattice, which is = A^(A^ -|- l)/2, 
iV^ A^(3A^ - l)/2, ~ ?,N + I, N"^ - {N - l)/2 for the triangle, rhombus, trapezoid, 
hexagon and rectangle shaped system, respectively. One can alternatively define the free 



energy density according to the actual geometrical area, which is different from the present 
definition by a trivial constant. 

With the BP algorithm, we obtain the free energy density / directly. The internal energy 
and specific heat are calculated by using a differentiation method 

+ A/3) -/(/?,- A/3) 



u 



2A/3 
I 

(A/3)2 



^ ^ ^2 /(/3c + A/3) + /(/3,-A/3)-2/(/3,) ^ 



In our calculation, A/3 = 10~^ is used. The analysis of error in the calculations of u and c 
has been discussed in 27|. The final estimation of the accuracy of the free energy, internal 
energy and specific heat is 10~^^, 10~^^, 10~^, respectively. 

The calculations were carried out for 105, 104, 103, 85, 127 systems, with linear size N 
varying from 30 to 2000, for the triangle, rhombus, trapezoid, hexagon and rectangle shaped 
triangular lattice, respectively. 

III. RESULTS 

A. Critical free energy density 

By fitting the finite size data, we find that the exact expansion of the critical free energy 
can be written in the following form, with k from 3 to 12 for the triangle and rectangle, k 
from 3 to 9 for the other shapes: 

p{N) fcornlnN + f2 ^ fk + k^nN 
J Joo + J surf g g / , ' \ / 

A;=3 

where p{N) is the perimeter, which equals to 3A^, AN, 5N — 1, 6N, (2 + \/3)N for the triangle, 
rhombus, trapezoid, hexagon and rectangle, respectively. This ex pan sion is different from 
that for the triangular lattice with periodic boundary conditions 18|, in which there's no 
surface, corner terms, logarithmic corrections, and only even k presents. 

The fitting method is the standard Levenberg-Marquardt method for nonlinear fit. The 



standard deviation (SD) is defined by a = yJ2iifi " tV^^'^Y I ('^rf ~ ''^z) "with fi the numerical 
data, /y**"* the value given by the fitting formula, the number of data used and nj the 
number of fitting parameters. For all cases, a reaches 10~^^. The accuracy is seen from the 
fitted bulk free energy density foo'- the worst fit among the five cases is for the rhombus. 



which yields foo = 0.8795853861615715170938962(7), and the best one is for the triangle 
yielding f^o = 0.87958538616157151709389605(3). These results coincide with the exact 



value /oo = 0.8795853861615715170938960283- • • (l7| in more than 24 decimal numbers. 



TABLE I. The fitted edge, corner free energy and /2 in Eq. 



shape fsurf fcorn /2 

triangle -0.1030776388340906553343(2) 0.1666666666666666667(4) 0.006804832446685952(4) 
rhombus -0.103077638834090655334(2) 0.14583333333333333(2) 0.1839728334687587(1) 
trapezoid -0.103077638834090655335(1) 0.14583333333333335(2) 0.2213952601428039(1) 
hexagon -0.103077638834090655335(2) 0.10416666666666668(3) 0.4758271527774812(2) 
rectangle -0.109078407337305392239(2) 0.125000000000000001(2) 0.22549835683994714(1) 

According to finite-size scaling, the surface correction term fsurfP{N) / S stems from free 
edges. As shown in Fig. [U^e), there are three bond directions for the triangular lattice. 
For the triangle, rhombus, trapezoid, and hexagon shapes, the edges are all along one bond 
direction. The edge free energy per unit length on these edges should be equal. For the 
rectangle shape, there are two edges along one bond direction and two other edges perpen- 
dicular to that bond direction in a zigzag way. Therefore the edge free energy per unit length 
along the two different directions can be different in principle. To be clear, we denote the 
surface free energy per unit length (the straight length, not the total length of the zigzag 
line) along and perpendicular to the bond direction as and f^rj, respectively. The 
surface free energy is thus fsurfP{N) = fgurfPi-^) ^'^^ triangle, rhombus, trapezoid, and 
hexagon shape, respectively. For the rectangle shape, it is fsurfP{N) = fsurfC^ + y/S)N, 
with fsurf = (2/|'„r/ + V^/snr/)/(2 + V^) , which is the mean surface free energy per unit 
perimeter, considering the straight length of the edges along the perpendicular direction is 
\/3N/2, which is used to define the perimeter. The fitted values of fsurf is given in Tab. [H 
For the rectangle case, we obtain 

/J^^^ = -0.103077638834090655334(2), f^^^.j = -0.116007497958656704304(2). (6) 

Cardy and Peschel showed in 3j that the presence of a corner of interior angle 7 along 
boundaries of typical size give rise to a logarithmic correction to the critical free energy 
density 

^(7) = _ (1)2)^]^ (J) 

Jcorn g 247r^ S ' 



where c is the conformal anomaly. The total corner free energy is 



fa 



q /(t/S) 
o J corn 



0.16666666. 



2(/, 
6/, 



(^/3) 
corn 



J cor' 



{27r/3) 



^ = 0.145833333... 



(277/3) 

corn 



^ = 0.104166666. 



J corn 



0.125 



for triangle 

for rhombus and trapezoid 
for hexagon 
for rectangle 



The fitted corner free energy fcom for the five shapes are listed in Tab. HI from which one 
can see that our results reproduces the CFT result very accurately. 

TABLE II. The fitted parameters for order /c > 3 in Eq. ([5]) for the critical free energy of the 
triangle and rectangle shaped triangular lattices. There is no logarithmic corrections except for 
the corner term, i.e., all = 0. 



shape triangle 



rectangle 



/3 
U 

h 
k 
h 
h 

/9 
/lO 

hi 

fl2 



0.117851130197757931(2) -0.0144959340650645(9) 





-0.00245523187923(4) 

0.001247090597(1) 

-0.0016255469(3) 

0.00124715(1) 

-0.0005633(2) 

-0.000553(3) 

0.00159(2) 

-0.050(4) 



0.0180683705511(1) 

-0.01237101066(1) 

0.0168641489(1) 

0.00302743(4) 

0.001335(1) 

-0.03133(3) 

0.0701(4) 

-0.161(2) 



In the previous study on the finite square lattice in a rectangle shape [27[, we estimated 
the corner free energy fcom = 0.125 ± 2.0^^° for the rectangle shape with various aspect 
ratios. Here we find the same result fcom = 0.125 ± 2.0 x 10^^*^ on the rectangle-shaped 
triangular lattice. This indicates that the corner term of the free energy is independent of 
the microscopic properties of the lattice. Therefore we proved the CFT prediction that the 
corner free energy is universal . However, our calculations are for the aspect ratio fixed as 
p = 2/\/3 at present work. Thus it is not enough to further verify Kleban's CFT predictions 



on the effect of the aspect ratio of the rectangle-shaped lattice 



2l|. 



TABLE III. The fitted parameters of Eq. ([5]) for the critical free energy. The third and fourth 
order logarithmic correction 1^,14 are zero. 



shape rhombus 


trapezoid 


hexagon 


/3 


DzL^^fil 11111111 9(^9^ 


-0.00210692865756(2) 


-0.06014065304055(7) 


fd 


-0.00810185186(2) 


01 363281 706l'2') 


0347222221 7l'7') 


J 


0.00381144(6) 


-0.00686168(6) 


-0.0180976(4) 


ffi 


-0.00215(1) 


-0.00894(1)0 


0.0150(1) 


f? 

J 1 


-0.0008(5) 


0.0316(6) 


-0.018(7) 


J 


—0 026f4) 


-0.113(5) 


-0.56(9) 


/9 


0.032(2) 


0.097(3) 


0.21(3) 


^5 


0.00346339(1) 


0.00494874(1) 


0.00999797(6) 




-0.003461(2) 


-0.008767(3) 


-0.00863(2)0 


h 


0.0024(1) 


0.0117(2) 


0.011(2)0 


h 


0.006(2) 


0.011(3) 


0.17(5) 


k 


0.007(5) 


0.070(8) 


0.9(2) 



For higher order terms, we give fitted coefficients in Tab. [TTl and lllli For the critical 
free energy of the triangle and rectangle shaped triangular lattices, there is no logarithmic 
corrections except for the corner term, i.e., all Ik = 0. For the triangle shape, the coefficient 
fd is determined to be zero since it is extremely small in our fits if this term is present, 
and the fitting result changes little if we discard this term. In contrast to the triangle and 
rectangle shape, for the other three shaped triangular lattices, we find higher order {k > 5) 
logarithmic corrections besides the corner term. Although these logarithmic corrections are 
very weak, with very small coefficients, our high accurate data indicate their existence. 



B. Critical internal energy density 

We fit the data on the critical internal energy with the following formula 

p{N)\nN , IniV ^j^ 

U Wjx) + Ugurf g ~v Ucorn ^ i / ^ Qk/2 ' ^ ' 

fc=l 

where p{N) is again the perimeter. In our fits, k is truncated to 4. Again, this formula 



is different from that for the triangular lattice with periodic boundary conditions 



18|, in 



which there's no logarithmic sur: 
bulk value Mqo is known to be 2 



"ace term, no corner term, and only odd k presents. The 



ITJ. Our fit of Moo is 2.0 ± 1.0 x lO^^^ for the five shapes. 



The other fitted parameter are given in Tab. IIVI 



TABLE IV. The fitted parameters of Eq. ([8]) for the critical internal energy per spin. 



shape triangle 


rhombus 


trapezoid 


hexagon 


rectangle 


Usurf -0.55132891(2) 


-0.551328867(8) 


-0.55132888(1) 


-0.55132891(2) 


-0.5513290(7) 


Ucorn —1.653995(3) 


-0.735055(7) 


-0.73506(1) 


1.10268(2) 


0.147707(6) 


ui 0.1358131(3) 


-0.7286080(2) 


-1.4553110(3) 


-3.4498392(7) 


-1.1246224(2) 


U2 -1.42006(1) 


-1.09359(3) 


-0.55400(4) 


1.71014(4) 


-0.95485(2) 


Us -0.70291(8) 


-0.2154(3) 


0.0302(4) 


-0.372(1) 


0.1049(2) 


U4 -0.0243(4) 


0.113(2) 


0.044(4) 


0.49(1) 


-0.030(2) 



The leading correction is due to the edges (or surface), which is unknown to our knowl- 
edge. We denote the surface internal energy per unit length of the edge along one bond 
direction by ul^^j, and that perpendicular to that direction by For the triangle, 

rhombus, trapezoid, and hexagon shapes, we have Ugurf = ""Lr/- -^^r the rectangle we have 
(2«Lr/ + ^^iur/)/ (2 + V3), and 

«L/ = -0.55132889(2), uj^^^ = -0.55132895(7). (9) 

This is an interesting result because it means that the surface internal energy per unit 
length is symmetric for an edge along or perpendicular to an arbitrary bond direction. 
Comparing this result with Usurf for the rectangle shaped square lattice in our previous 
workj^ll, where Usurf = 0.6366198(1), we find that the ratio is 1.1547005, which is very 
close^ 2/\/3 = 1.154700538 • ■ ■ . Because the exact value of Usurf for the square lattice is 
2/7r [11], we conjecture that the exact value of the surface internal energy for the triangular 
lattice is given by 

i/ = ujurf = -Vs/n. (10) 

Following the convention for the critical free energy, we write the coefficient of (In N) / S 
as Ucorn- We denote the corner correction by Umln, where 7 is the angle of the corner. Under 

10 



the assumption that Ucom is the sum of the comer's contributions, we have 



o (-/3) 
OLicorn 

I Uicorn \ Uicorn 



Q^(2-/3) 
\j LI corn 



Au 



(^/2) 



for triangle 

for rhombus and trapezoid 
for hexagon 
for rectangle 



From the Tab. llVt we obtain the corner term for the three angles 



= -0.551332(1), 



u 



W2) = 0.036927(2), 



u 



(2„;:/3) = 0.183781(3). 



In the previous study on the square lattice |27|, we obtained the corner term Ucom = 
—0.4502(1) for the rectangle with various aspect ratios. It is different from the present 
result Ucom = 0.147707(6) for the rectangle-shaped triangular lattice, which indicates that 
the corner term of the internal energy depends on the microscopic structure of the lattice, 
thus is not universal. 

The other parameters Bi, B2, B^, B4 are also estimated and listed in Tab. IIVI We have 
tried other forms of formula to fit the critical internal energy. The terms \n S/S^^"^, In S/S'^ 
are excluded considering the coefficients are extremely small. Moreover the standard devia- 
tions of the fits with these terms included are much larger than those without them. 



C. Critical specific heat 

The data of the critical specific heat are fitted using the following formula 

C= AolnN + Co + Csurf ^ h Ccorn— 7^ + 2^ "7;^, (12) 

k=l 

where p{N) is the perimeter, k is from 1 to 4 in our fits. Compared with the expansion for 



18| , there are additional logarithmic 



which reads An 



the triangular lattice with periodic boundary conditions 
surface term and corner term. 

The leading term A^lnN is known from the exact result 
^(ln3)2 ^ 0.4990693780 ■■■ . Our fit yields Aq ^ 0.499069374(5). The other fitted 
parameters are listed in Tab. IVl 

The leading correction p{N) In N/ S is caused by the edges. We denote the surface specific 
heat per unit length of the edge along one bond direction by cl^^^, and that perpendicular 

11 



to the direction by cj-^^f- triangle, rhombus, trapezoid, and hexagon shape, we have 

Csurf = cl^rf^ ^^le rectangle, we set Csurf = (2c^^^^ + V3cj^^j)/{2 + a/3), and find 



c 



surf 



0.166354(2), 



C" 



surf 



0.105462(2). 



(13) 



Note that this term is absent in the torus case [8] and not mentioned in the loiig strip case 



but exists in the cylinder case with Brascamp-Kunz boundary conditions 14j. Il5|. 



TABLE V. The fitted parameters of Eq. (I12p for the critical specific heat of the five shapes. 



shape triangle 


rhombus 


trapezoid 


hexagon 


rectangle 


Co -0.80424237(1) 


-0.60510331(1) 


-0.51995011(1) 


-0.2673395(1) 


-0.573388895(9) 


c,urf 0.1663558(3) 


0.1663531(5) 


0.1663520(5) 


0.166338(9) 


0.1380941(8) 


ci -0.187862(5) 


-0.036692(3) 


-0.10391(2) 


-0.2749(2) 


-0.08529(9) 


Ccorn 0.7484(1) 


0.3864(2) 


0.3857(3) 


-0.343(6) 


-0.1510(2) 


C2 0.6410(4) 


0.6443(8) 


0.615(1) 


0.58(2) 


0.6833(7) 


C3 0.339(2) 


0.177(6) 


0.04(1) 


-0.09(3) 


-0.057(5) 


C4 0.024(8) 


-0.12(4) 


-0.00(7) 


1.8(3) 


0.17(3) 



Following the convention for the critical free energy, we write the coefficient of (In N) / S 
as Ccorn- Wc dcuote the corner correction by Cmln where 7 is the angle of the corner. Again, 
under the assumption that the total correction is the sum of the corners, we have 



3 c, 



(27r/3) 
corn 



for triangle 

for rhombus and trapezoid 
for hexagon 
for rectangle 

From Tab. IIVI we obtain the corner contribution of the three angles. 



2 ("c^"/^) + c 

^ I ^corn 1^1 
'J^^corn 
^^^corn 



.(^/3) 



0.24948(3), 



.(^/2) 



-0.03775(5), 



c 



(27r/3) 



-0.057(1). (14) 



In the previous study on the square lattice [27|, we obtained the corner term Ccom = 
0.368(1) for the rectangle with various aspect ratios. It is different from the present result 
Ccorn = —0.1510(2) for the rectangle-shaped triangular lattice, which indicates that the 
corner term of the specific heat depends on the microscopic structure of the lattice, thus is 
not universal. 
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We have also tried other forms of fitting formula to fit the critical specific heat data. For 
example, we added the terms (In S)^/S, (In Sy/S in the fitting formula and found that their 
coefficients are extremely small. 

IV. CONCLUSION 

. Using the BP algorithm, we have studied the 2D critical Ising model on a triangular 
lattice with free boundaries. For five shapes, triangle, rhombus, trapezoid, hexagon and 
rectangle, the critical free energy, internal energy and specific heat have been calculated. 
We have proved the conformal field theory prediction of the corner free energy and have 
shown that the corner free energy, which is proportional to the central charge c, is indeed 
universal. For the edges parallel or perpendicular to the bond direction, the logarithmic 
edge corrections in the internal energy have been found to be almost identical, while these 
corrections in the free energy and in the specific heat have been found to be different. 
Comparing with the previous result on the square lattice in the rectangle shape, we have 
found that the corner internal energy Ucom and the corner specific heat Ccom for the rectangle 
shape arc not universal, i.e., the coefficient in front of In N/ S in the expansion of the internal 
energy and the specific heat are different for the square lattice and the triangular lattice. 

We have also found that there exist logarithmic corrections in higher orders, say, there 
are terms In N/S^^'^, InN/S^, - • ■ in the critical free energy for the rhombus, trapezoid, and 
hexagon shapes. However these terms are absent for the triangle and rectangle shapes. This 
should be an interesting subject to be further investigated. 
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